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Abstract. A nano-system in which electrons interact and in contact with Fermi leads gives rise to an 
effective one-body scattering which depends on the presence of other scatterers in the attached leads. 
This non local effect is a pure many-body effect that one neglects when one takes non interacting models 
for describing quantum transport. This enhances the non-local character of the quantum conductance 
by exchange interactions of a type similar to the RKKY-interaction between local magnetic moments. 
A theoretical study of this effect is given assuming the Hartree-Fock approximation for spinless fermions 
of Fermi momentum fcp in an infinite chain embedding two scatterers separated by a segment of length 
L c . The fermions interact only inside the two scatterers. The dependence of one scatterer onto the other 
exhibits oscillations of period 7r/fcF which decay as 1/L C and which are suppressed when L c exceeds the 
thermal length Lt- The analytical results given by the Hartree-Fock approximation are compared with 
exact numerical results obtained with the embedding method and the DMRG algorithm. 

PACS. 71.27.+a Strongly correlated electron systems; heavy fermions - 72.10.-d Theory of electronic 
transport; scattering mechanisms - 73.23.-b Electronic transport in mesoscopic systems 



1 Introduction 

The coupling of nano-objects via conduction electrons was 
discovered long ago, in the case of spins of magnetic ions, 
or of nuclei, which results indirectly from the interaction 
of such spins with those of conduction electrons in met- 
als. After tracing out the degrees of freedom of the con- 
duction electrons, one gets an effective spin Hamiltonian 
characterized by an oscillatory long range interaction, the 
RKKY-interaction [TJdlGIlH] > which plays a crucial role in 
understanding the large variety of possible ordered spin 
structures in magnetic crystals. If Jj is a local magnetic 
moment in a metal, the conduction electrons give rise to 
an interaction energy between these moments, which can 
be described by an Hamiltonian of the form: 



H 



RKKY 



(1) 



The coupling term between two moments separated by a 
distance Ra behaves as 



Jij OC 



2kpRij cos(2fcF-Ry) — sin(2/cFi?ij) 



(2) 



in d = 3 dimensions, being the Fermi momentum of 
the conduction electrons. In d = 1 dimension, this gives a 



long range 1 / R interaction with oscillations of periodicity 

7r//CF- 

We show in this work that a similar phenomenon char- 
acterizes also the quantum conductance of nano-systems 
in which electrons interact, coupled via metallic wires. 
This phenomenon is very general and does not require 
to include the spin degrees of freedom. Combining Lan- 
dauer's formulation [5] of quantum transport and the 
Hartree-Fock approximation, as reviewed in Ref. [6l[7j , we 
will show that the scattering matrix of an interacting sys- 
tem depends on what is embedded at a distance R of the 
interacting system in the attached leads, this dependence 
decaying as 1/R d with oscillations of periodicity 7r/fc F . 
This non local character of the quantum conductance is 
another example of the effect of indirect exchange inter- 
actions between interacting nano-systems via conduction 
electrons, as the RKKY-interaction between local mag- 
netic moments. 

To study this phenomenon, we take an infinite chain 
where spinless fermions do not interact outside two iden- 
tical regions where the scattering is a pure many body 
effect due to Coulomb repulsions. The chain is described 
by a one dimensional tight binding model with nearest 
neighbor hopping th — 1. A nearest neighbor repulsion 
of strength U acts only between two consecutive sites, in 
two regions which are connected by L c sites where the par- 
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tides do not interact. This model is simple enough to be 
analytically solved at a mean field level, using the Hartree- 
Fock (HF) approximation. To simplify the calculations, we 
add positive compensating potentials which exactly can- 
cels the Hartree terms of the HF equation. But the ex- 
change term modifies the hopping term coupling the two 
internal sites of each scatterers, in such a way that it takes 
a value v instead of th- In the HF approximation, v is given 
by a self-consistent equation. Having v, it is straightfor- 
ward to obtain the scattering matrix Si (v) of one scatterer 
at the Fermi energy of the infinite non interacting chain. 
Using the combination law of the one-body 5-matrices in 
series, one can get the scatteringmatrix ^(w), and hence 
the dimcnsionlcss conductance g 2 F (in units of e 2 /h) of 
these two many-body scatterers in series, in the HF ap- 
proximation. 

For such a system where two identical many-body scat- 
terers are coupled by a perfect wire of L c sites where the 
electrons do not interact, our main result is to show that 
the value of v characterizing one scatterer differs from its 
value when there is no second scatterer by an oscillatory 
term of period 7r/fcp which decays as l/Lc. In certain lim- 
its, v can be given in a simple form as a function of U, L c 
and kp. For instance, for a half- filled chain (fcp = 7r /2) 
and weak interactions, the effective hopping term v reads 

2 (-i) l - / i y 

V ~ + 2-tt + 2tt/U + L c \2-tt + 2tt/UJ + "' 

In this limit, this expression shows how the effective hop- 
ping term v characterizing a single scatterer decays as 
1/L C towards its value when it is not in series with an- 
other one, with the even-odd oscillations characteristic of 
a half-filled chain. Using Eq. ©, it is straightforward to 
show how the scattering matrix Si(v) of a single scatterer 
is modified by the presence of a second scatterer. We un- 
derline that this non local effect is a pure many-body effect 
that one neglects when one takes non interacting models 
for describing quantum transport. This non local effect 
was first numerically discovered in a previous work J8| , us- 
ing the embedding method piT0l[TT|[T^[T^[T4l[l5l[T6] and 
the DMRG algorithm [TTlfTS] valid for one dimensional 
fermions. In this work, we give a simple theory of this ef- 
fect based on the HF approximation, which turns out to 
qualitatively describe this non local effect for all values of 
U, including its suppression in the limit when U — » oo. 
The HF approximation becomes quantitatively accurate 
for small strengths U of the interaction. Moreover, we 
will also show that this non local effect vanishes when the 
length L c of the coupling wire exceeds the thermal length 
I/t characterizing free fermions in one dimension. 

The paper is organized as follows: 

In section [21 we consider a single scatterer with a near- 
est neighbor repulsion of strength U embedded in an in- 
finite chain. In the first sub-section, the HF equation of 
this simple model is written, leading us to study a chain 
where the hopping term between the two central sites is 
equal to v instead of th = 1. This one body model is 
solved in the second sub-section, allowing us to obtain the 
implicit equation giving v in the HF approximation. The 



system being symmetric upon reflection, we use scattering 
phase shifts and Friedel sum rule for this purpose. In the 
third sub-section, the conductance of this single scatterer 
is studied as a function of the strength U of the nearest 
neighbor repulsion, the HF behavior being compared to 
the exact results given by the embedding method and the 
DMRG algorithm. In the fourth sub-section, the Friedel 
oscillations of the particle density around the scatterer are 
described. In a last sub-section, a correlation function in- 
side the attached leads is calculated at a distance p from 
the scatterer, which will be useful for describing the case 
of two scatterers in series. This function is shown to decay 
as 1/p with oscillations of periodicity n/k-p towards an 
asymptotic value characterizing the chain without scat- 
terer. 

In section [3l we study the conductance of two scatter- 
ers in series, coupled by a scattering free wire of L c sites. 
In the first sub-section, simple analytical expressions are 
given in the limit where the strength U of the interaction 
acting inside the two scatterers remains small. Notably, we 
show that the effective hopping term v characterizing each 
scatterer differs from the value obtained in the section [5] 
by a correction which decays as 1/L C with oscillations of 
periodicity n/k-p. In the second sub-section, the Hartree- 
Fock equation is solved exactly for arbitrary values of U, 
allowing us to show that the weak [/-expansion assumed 
in the first sub-section remains valid till U w th. The effec- 
tive hopping term v characterizing a single scatterer being 
modified when it is in series with another, the implication 
of this non local effect upon the conductance is illustrated 
in the third sub-section, the HF curves describing the con- 
ductance oscillations of the two scatterers in series being 
similar to the curves numerically obtained in Ref. [5]. 

In section |4l the results for two scatterers in series 
given by the HF theory are compared to the exact numer- 
ical values given in Ref. [8|. The HF theory turns out to 
give a good qualitative description of the non local effect, 
which becomes quantitatively accurate when U < th for 
the very small scattering regions which we have consid- 
ered. 

In section 03 we show that this non local effect is sup- 
pressed at a temperature T, when the length L c of the 
coupling wire exceeds the thermal length Lx oc vp/T of 
free fermions in one dimension, vp being the Fermi veloc- 
ity. 

We conclude in section [6] by a summary of the main 
results, underlining their relevance for a theory of the non 
local character of quantum transport measurements, and 
suggesting straightforward extensions of this theory out- 
side one dimension. 



2 Transmission through a single many-body 
scatterer 

2.1 Microscopic model and exchange energy 

To study the indirect exchange interaction via conduction 
electrons between nano-systems in which the electrons in- 
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teract, we begin to study the simplest many-body scat- 
terer, taking a tight-binding model of spinless fermions 
which do not interact, unless they occupy the two central 
sites and 1 of an infinite chain, which costs a nearest 
neighbor repulsion energy U. Assuming the Hartree-Fock 
approximation, this leads us to study an analytically solv- 
able one body model where the hopping term between 
the two central sites is modified by the interaction, mod- 
ification which has to be calculated self-consistently. The 
many-body scattering system is described by an Hamilto- 
nian H — Hki n + Hint ■ The kinetic and interaction terms 
respectively read 



■Mr) -1 



Hkin— ^h(CpC p _ 1 + C p _ 1 C p ) 

p— — oo 

H lnt = U[m- V+] [n - V+] . 



(4) 



The hopping amplitude fh = 1 between nearest neighbor 
sites sets the energy scale, c p (cD is the annihilation (cre- 
ation) operator at site p, and n p = c p c p . The potential V+ 
is due to a positive background charge which exactly can- 
cels the Hartree term of the HF equation. The conduction 
band corresponds to energies —2 < E = — 2cosfc < 2 (fc 
real). H is invariant under reflections (p— 1/2 — > — p+1/2) 
and exhibits particle-hole symmetry if the chain is half 
filled. In this case, V+ = 1/2, the Fermi momentum 
k-p = 7r/2 and one has a uniform density without Friedel 
oscillations around the central region where the fermions 
interact. 

In the HF approximation, one assumes a variational 
ground state which is a Slater determinant of one-body 
wave- functions ip a (p) of energies E a < Ep. Since in our 
model the negative charge inside the scatterer is exactly 
compensated by a positive background charge, the Hartree 
term is cancelled and we have just to take into account the 
exchange term in the Hartree-Fock equation [T5|I^D] giving 
the tpaip): 

-V Q (p + l)-^(p-l)-^0^(P') = ^a(p). (5) 

P' 



The exchange term 

E,3<E F 

with the taken nearest neighbor repulsion 

U p , p i = U {8 p fi8 p r,i + S p ,iS p /,o} , 



(6) 



(7) 



is very simple. It is zero, excepted between the two cen- 
tral sites where the fermions interact, where it yields an 
increase of the strength of the hopping term coupling the 
two central sites by an amount 

tZ(U) = U J2 rM)^(0) = u(c{c o ). (8) 

E a <E F 

The Hartree-Fock equation describes a tight-binding 
model which is represented in FigdJ where the hopping 



) — 







— • — ( 

i 1 



Fig. 1. Effective one body model obtained assuming the 
Hartree-Fock approximation for a single many body scatterer 
(Hamiltonian @). The hopping term (indicated above) v(U) 
between the sites at p = and p = 1 (indicated below) depends 
on U and fcp. 



term between the two central sites is no longer equal to 
th = 1, but takes an interaction dependent value v, which 
is given by an implicit equation: 



HF/ 



(9) 



for t h = 1. 



2.2 Scattering phase shifts and density of states 

The effective one-body model being symmetric upon the 
reflection (p— 1/2 — > —p + 1/2), its Hamiltonian has even 
and odd standing- wave solutions V'fe(p) an d ipk(p)i which 
can be written [21j inside the conduction band as: 



(10) 



^(p)=y|cos (k(p-±)-6a(k) 
4>l{p)=J\s^(k{p- 1 -)-5 l {k) 



2 / I 
-cos I k(p- -) + S (k) 



(11) 



at the left side of the scatterer (p < 0) and 

1 2 / 1 
^(p)=y-jrsin U(p- -)+Si (fc) 



at its right side (p > 1). The normalization factor y jr 
corresponds to a chain of length L — * oo. The scatter- 
ing when u / 1 gives rise to two phase shifts 6o(k) and 
<5i(fc). Writing the Schrodinger equation inside the scat- 
terer (sites and 1) for the even and odd solutions, one 
gets: 



-(2cos/fc)V'°' 1 (0) 
-(2cosfc)V^' 1 (l) 



^(o)-^), 



(12) 



which yields the following expressions for the even and 
odd phase shifts: 



and 



. ,,, v — 1 
tan d (k) = — — rcot 

v + 1 V 2 



r . 1 — v ( k 

tanoi Ik) = tan — 

n ; w+1 \2 



(13) 



(14) 



4 



Y. Asada et at: Indirect exchange interaction 



In addition, a value of v > 1 gives rise to two bound 
states located in the central region with energies outside 
the conduction band. The first one of energy 



E hs i = -(v + v' 1 ) 
has an even wave-function given by 



,-b-l/2| 



V-bsib) = N — — 

while the second of energy 

£b S 2 = v + ir 1 
has an odd wave-function given by 



(15) 



(16) 



(17) 



(18) 




0.75 



0.5 - 



0.25 



U 

Fig. 2. Conductance of a single scatterer as a function of 
the interaction strength U for fcp = n/2. The solid line gives 
the HF behavior (Eq. (J25J ) . The circles are the exact results 
obtained with embedding method and the DMRG algorithm. 



When a scatterer is introduced in the chain (v ^ 1), 
this yields a correction Sp(E) to the density of states p(E), 
which is given [22 for the even and odd components inside 
the conduction band by: 



1 dSoM 

7T dE 



8p(E) satisfies [22] Friedel sum rule: 
S (E F ) + 5i(E F ) 



Sp{E)dE = vL s 



(19) 



(20) 



for a scatterer of length L$ embedded in a chain with a 
uniform filling factor v. This implies that 



Sf'Sp(E)dE Ls 



J^p(E,v = l)dE L 



(21) 



when L — > oo: The change Sp(k) of the density of real mo- 
menta k (inside the conduction band) vanishes in the limit 
of infinite lead length. When L — > oo, the only change in 
the density p(k) due to the scatterer comes from the bound 
states -0bs of imaginary momenta K — ik which are oc- 
cupied outside the conduction band at zero temperature. 
For an arbitrary function F{k) at zero temperature, this 
gives the relation: 

fc<fc F JV bs 

the last term being a sum over the occupied bound states. 



2.3 (c}c ) and conductance gi F (v) 

To obtain the implicit equation giving v, we need to calcu- 
late (ejeo), which is the sum of contributions due to 
the conduction band and Alf due to the occupied bound 



states. Assuming that Ep < 2, only the even bound state 
is occupied and: 



c{ c 



A 1 ' - 




C(i)#(o)| 

Av cos k sin 2 k 
1 + v i - 2v 2 cos(2fc) 

2v sin fcp \ sin k-p 



(23) 



• arctan 



4 1,0 
^Sl 



2tt 

= ^ sl (l)Vbsl(0) 



- 1 



1 



Using this, one can calculate v as a function of the 
interaction strength U and of the Fermi momentum kp by 
solving the implicit equation 



v = 1 + U(c\c (v)). 



(24) 



Once the change of the effective hopping term v between 
the two central sites is obtained, it is straightforward to 
determine the transmission coefficient ti(v). At zero tem- 
perature, the Landauer conductance ff 1 F (w) (in units of 
e 2 /h) of the central region where the electrons interact is 
given by this effective one-body transmission coefficient 
|^ F ( W )| 2 - Using the Landauer formula, one gets the trans- 
mission coefficient and the dimensionless conductance 



v(e 



2i/vi 



1) 



v*- — c 

^ F («)l 2 



2iA F 



4w 2 sin 2 kp 



(25) 



v 4 - 2v 2 cos(2fc F ) + 1 ' 



The behavior of as a function of the interaction 
strength U, obtained assuming the Hartree-Fock approx- 
imation (Eq. (|25|) ) is shown in Fig. [2] for fcp = tt/2. An 
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A 0.04 




P 

Fig. 3. Particle density n p = (c P c p ) at sites p around a scat- 
terer with an interaction strength U = 0.4. For &f = n/32 
(lower part of the figure): Contributions of the bound state 
(dashed line) and of the conduction band (dotted-dashed line) 
to rip and Friedel oscillations of n p around the average fill- 
ing 1/32 (solid line). For &f = tt/2 (upper part of the figure): 
Contribution of the conduction band (dotted-dashed line) and 
total uniform density n p — 1/2. 



accurate value for g 1 ^ can be obtained using the embed- 
ding method and the DMRG algorithm, as introduced in 
Refs. [9lfT0], Using this exact method, we have also cal- 
culated <7 X F at different interaction strengths U. The data 
are presented in Fig. [21 showing that the HF approxima- 
tion is a good approximation for a very short interacting 
region of moderate interaction strength U. 



2.4 Friedel oscillations of the density outside half-filling 

For a half-filled chain, the system has particle-hole sym- 
metry, and the density n v = (cpc p ) must be equal to 1/2 at 
each site p. This is due to two opposite effects which com- 
pensate each other: a decrease of the contribution of the 
conduction band to the density that one can expect when 
there is a local repulsion acting inside the scatterer, and 
the contribution of the interaction-induced bound state to 
the density. These contributions are plotted separately in 
Fig- El For half-filling (upper part of Fig. [3]), one gets an 
exact compensation of the two opposite effects for having 
a uniform density, and particle-hole symmetry is satisfied. 
Outside half- filling (lower part of Fig. [3]), one gets large 
Friedel oscillations of the density around the scatterer, as 
shown in Fig. [3] for fep = tt/32 (filling factor 1/32). Since 
the Hartree term is exactly compensated in our model, 
those oscillations are only due to the exchange energy. 



2.5 



p+l'-Ti 



outside the scatterer 



The conductance g 1 F (v) was obtained assuming infinite 
perfect leads outside the central scattering region. To 
know to what extend this condition does really matter, 



0.34 
0.33 

-J 0.32 

+ 

J? 0.31 
0.3 
0.29 



■ » o o 

o ° o ° > ° o ° o 
o 



10 



15 



20 



25 



30 



Fig. 4. Oscillatory decay of (c P+1 c p ) towards its asymptotic 
value 1/tt as a function of p for /cf = n/2, U = 0.1 (circles) 
and U = 0.5 (diamonds) obtained by using ([28} -([30]). 



we calculate (c P+1 c p ) outside the two sites where the elec- 
trons interact. When p > 1, 



J r \ — A p+1 ' p 4- A p+1 ' p 
c p+l P / — cb ^ ^bsl ' 



(26) 



where the contributions A^ 1 ' 11 of the conduction band 
and A^^' p of the bound state to (Cp +1 c p ) read: 



dk 



— {cos k — G(v, k)} 



G(v,k) =(v 2 -l) 



v 2 cos(2fcp — k) — cos(2fcp + k) 



(27) 



1 + v A - 2v 2 cos(2fc) 



/ip+i.p 



V Z - 1 



-2p-l 



After integration, one obtains: 
sin &f v 2 — 1 



c p+i c p 



X(k F ,p,v) 



7T 7r(2j5+l) 

where the function X(kp,p,v) is defined as 
X(k^,p 1 v) = lmF(kF,p,v) 



(28) 



(29) 



F(k F ,p,v) ^F 1 {l,±+p,l+ P ye 2U *)e ik *<- a * +r >, 
2-Fi(a, 0, 7, z) being the Gauss hypergeometric function 

(a)„(/3) n 



2 F 1 (a,f3,j,z) = ^2 



(30) 



with (x) n = Yim=o( x + TO )- Using an expansion of the 
Gauss hypergeometric function valid in the limit p — > oo: 



(31) 
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o o-o o 



-o— 



■ + 1 



Fig. 5. Effective one body model describing two identical 
many-body scatterers (Hamiltonian (|33[0 coupled by L c sites 
where the electrons do not interact. 



study this effect in the limit when U — > 0, before solving 
exactly the Hartree-Fock equation. We give only the re- 
sults for L c even. The extension to the case where L c is 
odd is straightforward. Moreover, for fc F = ir/2 and odd 
values of L c , the conductance of the two scatterers in se- 

7T /2 

ries g 2 = 1, independently of v. Because of this, we just 
need v for the even value of L c , at half- filling. 



one obtains for the behavior of (c P+1 c p } far from the scat- 
terer a simpler expression: 



c p+i c i 



If.. , H^(P) 



(32) 



where H£ (p) is an oscillatory function given by 



hUp) 



{v 2 - 1) [-v 2 sin(k F (2p - 1)) + sin(fc F (2p + 1))] 
1 + v 4 -2v 2 cos(2fc F ) 



Outside the region where the electrons interact, 
(Cp +1 c p ) exhibits oscillations of periodicity ir/kp which 
have a slow power law decay towards an asymptotic value 
sinfcp/7r. Those oscillations are illustrated in Fig. [4] for 
fcp = 7r/2 and two values of U, 



3 Transmission through two many-body 
scatterers coupled via conduction electrons 

We now consider an infinite tight binding chain with the 
same kinetic Hamiltonian than defined in Eq. ([4]) , but with 
an interaction Hamiltonian: 



Hint =U(nt^ +2 - Vi+){til. +1 - V 2+ ) 
+U{n_r^ - V 2+ ){n_^_ 1 - V 1+ ). 



(33) 



This Hamiltonian describes two nano-systems where 
the electrons interact, as previously studied, which are 
coupled via an ideal lead of L c sites where the electrons 
do not interact. Two positive background potentials Vi+ 
and V 2 + are introduced to compensate the Hartree terms 
of the HF equation. The total Hamiltonian is symmetric 
under the reflection p — 1/2 — > — p + 1/2. At half-filling, 
one has particle-hole symmetry, V\+ = V 2 + = 1/2 and 
the density n p = 1/2 is uniform. As before, the exchange 
leads to a modified effective hopping term v in the HF- 
approximation. Because of reflection symmetry, this mod- 
ification must be the same for the two scatterers. The 
corresponding one-body model is sketched in Fig. [5] How- 
ever, when there are two scatterers in series, the value of 
v characterizing each scatterer becomes different from the 
value of v obtained in the previous section when there is a 
single scatterer. This is due to the indirect exchange inter- 
action which takes place between two scatterers coupled 
via conduction electrons. This indirect exchange interac- 
tion gives rise to an effect upon the quantum conductance 
of two nano-systems in series, which vanishes only if the 
length L c of the coupling wire becomes infinite. We first 



3.1 Expansion in the weak interaction limit 

For writing the self-consistent equation giving the effec- 
tive hopping term v of one scatterer in series with an- 
other, we need to calculate the ground state expectation 
value of (cl c CL a . . (v, v))2 inside one scatterer when an- 

other identical scatterer is located at the sites — kf- — 1 
and — 2 ■ This value depends of the two modified hop- 
ping terms which characterize each scatterer, and which 
are equal because of reflection symmetry. In the limit of a 
weak interaction strength U, v — ► 1 + (v > 1) and one can 
expand: 



(ck +2 %+iM) 2 = (<4 +2 c ¥+1 (m)) 2 

+ j|-C 2 (L c ,<;)} +0(( V -1) 2 ) 



(34) 



where 
C 2 (L c ,v) 



CV Cl c . 



1(1. «)) +(4 +2 ^ +1 M)X 



We note that the above expansion involves only terms with 
a single scatterer, for which one can write 



TT+2 2 + 



-^+2 ^ + 



^M^^e-tcoU)) 
1 («,!)) =(4c o (t0) 



(35) 



* +a c^ + i(l,«))_=(4„+3 c ie+a(w) 



Using Eqs. and (|2^]l. one gets when v — > 1 (v > 1) 



dv 



c\co(v) 



1 sin fcp 

2 V 



3 2XV *-^ 



dv 



(36) 



7T(2p+ 1) 



the function X(kp,p,v) being defined in Eq. (|29l) . 

In the weak interaction limit, the effective hopping 
term v characterizing each of the two scatterers in series 
is given by the self-consistent equation: 



l-v^-U \- + (v-l)^-C 2 (L c ,v) 

IT OV 



(37) 



)->!+ 



where C 2 {L c ,v) = (c\cq(v)) + (c f Lc+3 c Lc + 2 {v)). One even- 
tually obtains: 



, n(2-U) 2A(fc F ,£ c + 2,l) . 
v ps 1 + <, — — — h sin kf 



2U 



2L r 



(38) 
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when U — > 0. One can see that the indirect exchange in- 
teraction gives rise to a correction which decays with a 
power law as the length L c of the coupling wire increases. 
This is when L c — > oo that the scatterers become decou- 
pled, and characterized by a value for the effective hopping 
term which coincides with the value given by Eq. (|24| for 
a single scatterer in the limit U — > 0: 



v(L c = oo,U -> 0) w I + 



2sinfc F - 7r + 27r/[/' 



(39) 



Equation (|38p can be written in a simpler form for a 
half-filled chain. Expressing X(kp, L c + 2, 1) for fcp = tt/2 
one obtains for the effective hopping term the equation 
^ given in the introduction. 



3.2 Exact solution of the Hartree-Fock equation 

Let us solve now the HF-equation for v without assuming 
that v — > 1 + . The effective one-body model described in 
Fig.[5]has even and odd standing- wave solutions ip®(p) and 
%jj\{p), which can be written inside the conduction band 
as: 



$(p)=^/f cos (k(p-±)-5 (k^j 



at the left side of the two scatterers (p < — % — 1) and 



(40) 



$(P) 



2 ( I 

-cos ( fc(p- -) +(5 (fc) 

|sinffc(p-±)^ A, 



(41) 



at the right side of the two scatterers (p > %-+2). Between 
the two scatterers (—4/ < P < 4r + 1), the even and odd 
standing- wave solutions ip®(p) and ip\{p) read 

ipkip) =\J sin ( k (P - \)) > 



^■aocos (fc(p- i) 



(42) 



The expressions for the factors a and a\ are given in 
Appendix. 

When v I, the scattering gives rise to two phase 
shifts <5g(fc) and §i(k), which are given by 



tan 5o(k) 



and 



tan Si (fc) 



respectively. 



(t; 2 - 1) (cos(/c(L c + 2)) + cos fc) 
(v 2 - 1) sin(fc(L c + 2)) + (v 2 + 1) sin A: 



(1 - v 2 ) (- cos(k(L c + 2)) + cos fc) 
(u 2 - 1) sm(k(L c + 2)) - (v 2 + 1) sin A: 



(43) 



(44) 



In addition, a value of v > 1 can give rise to four bound 
states located around the scatterers with energies outside 
the conduction band. We just write those below the con- 
duction band (Et, s < —2). There is an even bound state 
of energy i?bsO = — 2cosh.Koj Ko being the real solution 
of the equation: 

v 2 + v 2 exp(K (L c + 1)) = 1 + exp(^o(£ c + 3)). (45) 

Its wave function reads 

^bso(p) = A> exp \Ko(p + ^y^YJ (46) 
at the left side of the two scatterers (p < — 4/ — 1), 

^ bs0 (p) - A exp (-K (p - ^y^)) (47) 
at the right side of the two scatterers (p> 4/ + 2), and 



iphso{p) = A obo cosh ( K Q (p 



(48) 



between the two scatterers (— 4/ < p < 4/ + 1). The 
expression for the factor bo is given in Appendix. 
If L c is large enough, the equation: 

v 2 - v 2 exp(K 1 (L c + 1)) = l-exp(Kx(L c + 3)) (49) 

has a real solution for K\ . In this case, there is also an odd 
bound state below the conduction band of energy £?bsi = 
— 2coshJ*Ti. Its wave function is given by 

i> hs i(p) = -A 1 ex V (K 1 (p+^±±)) (50) 

at the left side of the two scatterers (p < — 4/ — 1), 

L c + 3, 



■0bsi(p) = M exp -Ki{p 



(51) 



at the right side of the two scatterers (p > 4/ + 2), an d 



ipbsiip) = Aibi sinh ( K\{p - -) 



(52) 



between the two scatterers (— 4f < p < 4/ + 1). The 
expression for the factor 6i is given in Appendix. 
The condition 



I^bs0,l(p)| 2 = 1 



(53) 



gives the normalization factors Aq>^. The obtained expres- 
sions are somewhat involved and given in Appendix. 
Using Eq. (T2"2"|) , one can calculate 



2 2tt 



(54) 
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Fig. 6. Two scatterers in series: Value v of the effective hop- 
ping term as a function of U. v characterizes each of the scat- 
terers for L c — 4 and = n/2. The solid line gives the exact 
HF value obtained from Eq. (154[) and the dashed line gives the 
approximated value (Eq. (|38[) ). 

1.0326 



1.03255 



s> 1.0325 



1.03245 



° 00 °°°°iioYoo"oooo 



OOOOOOOOOOOOOOOOOO 



1.0324 

20 40 60 80 100 

L c 

Fig. 7. Two scatterers in series: v as a function of even values 
of L c for U — 0.1 and kp = tt/2. The circles give the exact 
HF value obtained from Eq. (I54|l and the dashed line gives 
the approximated value (Eq. |38)| ). The even-odd oscillations 
characteristic of &f = n/2 are not shown, v for odd L c being 
not plotted. 




150 



Fig. 8. Two scatterers in series: Friedel oscillations of the den- 
sity n p = (cfpCp) for U — 0.4, L c — 100 and &f = n/32. 
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Fig. 9. Two scatterers in series: total conductance g 2 F (L c , U) 
as a function of L c for U = 0.4 and &f = n/32. Only the values 
given by the exact solution of the Hartree-Fock equation where 
L c is even are plotted. 



around the two scatterers. Outside half-filling, our model 
exhibits Friedel oscillations of the density which are illus- 
trated in Fig. E for U = 0.4, L c = 100 and k F = tt/32. 



for p = kj- + 1 to obtain the HF equation giving v for 
two scatterers in series. The above integrals have been 
calculated using Mathematica. 

We compare in Fig. [S] and Fig.[7]the exact HF value of 
v obtained with formula (|54|) to the approximated value 
given by formula (|38p in the limit v — * 1 + . One can see 
that the approximated value is accurate enough for the 
values of U where one can trust the HF approximation. 
Fig. [7] shows the effect of indirect exchange interaction 
upon the value v of the effective hopping, and how this 
effect disappears when L c — > oo, for even L c only. 

3.3 Density oscillations and quantum conductance for 
two many-body scatterers in series 

Once the self-consistent value for v is obtained, one can 
calculate the Friedel oscillations of the density n p = (c|cp) 



The transmission coefficient ij* of a single scatterer as 
a function of v and fcp is given by Eq. (p?5|) . r\ ¥ denoting its 
reflection coefficient, the transfer matrices M^ F and M^ F 
through a single scatterer and the coupling lead read: 

M k F _( (i/tin* r^/t<?\ (55) 

and 

From the matrix = M^ F M^M^ F , and current 

conservation (1 = |^ F | 2 + |^i F | 2 ), one obtains the trans- 
mission coefficient i 2 P an d the dimensionless conductance 
of two scatterers in series, coupled by L c sites where 
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0.913 
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2(H) 



Fig. 10. Two scatterers in series: Minimum values of 
52 B (Lc,U) for the successive conductance oscillations occur- 
ring when U — 0.4 and fcp = tt/16, underlining the power law 
decay of their amplitudes as a function of L c . The dashed line 
is a numerical fit 0.91324 - 0.05267/L c . 



the electrons do not interact: 

2iv 2 e 2ikF sin 2 kp 



d kF (v) 



4v 4 sin 4 kp 
" \d k *(v)\ 2 ' 



(57) 



where 



d kF {v) 



=e _ifeF sin[fc F (L c + 3)] - 2v 2 sm[k F (L c + 2)} 



-v 4 e ikF sm[k F (L c + l)}. 



The presence of L c dependent corrections to v shows 
that the transmission coefficient t\ F of one scatterer de- 
pends on the distance L c from the other scatterer. Ac- 
cordingly, this non-local effect affects the conductance g kF 
of two scatterers in series. We show in Fig. [9] the total 
conductance g 2 F of two scatterers in series for a low fill- 
ing factor (kp — it/32), obtained by the exact solution of 
the Hartree-Fock equation. One can see the large period 
n/k F = 32 of the conductance oscillations, though the 
values for odd L c are not plotted. The conductance oscil- 
lations are larger when L c is small, due to the effect of 
the exchange energy. The l/L c -decay of the conductance 
oscillations towards its asymptotic L c -independent value 
is underlined in Fig. [T0l obtained with kp — 7r/16. 



4 Comparison with exact DMRG results 

Exact values of the conductance of a one dimensional scat- 
terer in which electrons interact can be obtained using 
the embedding method and the DMRG algorithm, as ex- 
plained in previous works [8l[9T [TUl[TT) . We have compared 
in Fig. [5] the HF results and the exact DMRG results for 
a single very short scatterer. The difference was negligible 
for small values of U. Nevertheless, the HF values differ 
[7] more and more from the exact values when the size 



of the scatterer in which the electrons interact increases. 
The difference should become more pronounced for two 
scatterers in series. We study the ability of the HF ap- 
proximation to describe two interacting nano-systems in 
series for a half filled chain (kp — it/ 2). 

If the scattering matrix of one scatterer is not modified 
by the presence of other scatterer as for non-interacting 
systems, the conductance g^ 2 of two scatterers in series 
for kp — it/ 2 shows even-odd oscillation as a function of 



the size L c of the coupling wire: g 2 
and is given by [H] 



tt/2 



tt/2 
92 



1 when L r is odd, 



(58) 



when L c is even. Here g^ 2 is the conductance of a chain 
with a single scatterer. In the presence of electron-electron 
interaction, however, the scattering matrix of one scatterer 
can be affected by the other scatterer as shown in the 
previous subsection. 

In Ref. [8], the exact values of g^ 2 and g^ 2 were 

obtained separately, using the embedding method: g[^ 2 
being calculated for an infinite chain embedding a single 
scatterer, g^ 2 for an infinite chain embedding two scatter- 

77 /2 

ers. It was found that g 2 — 1 if L c is odd, as predicted. 

But if L c is even, g^ 2 was related to g^ 2 by formula (fS"5)) 
only in the limit L c — > oo. For small sizes L Cl formula 

7T/2 



(|58| , with g 1 obtained for a single scatterer surrounded 
by infinite leads without other scatterers, overestimates 
<?2^ 2 by an amount 

A{U,k F ) 



8gl' 2 (L c 



(59) 



characterized by a function A(U, kp) given in Fig. [Til for 
fc F = tt/2. 

Using the HF approximation, we have shown that 
the parameter v characterizing a single scatterer becomes 
modified if another scatterer is put in series. Hence, the 

7T /2 

conductance g 1 of one scatterer in series with another 
differs from its value when it is alone. If one ignores this 
difference, as previously using the DMRG algorithm, tak- 

77 /2 

ing for g x its value without the second scatterer and 

using formula (|58p . one overestimates also g^ 2 by an 
amount which is described by formula (|59p . but with a 
different function A(U, kp) given in Fig. QT] One can see 
that the HF approximation reproduces qualitatively the 
DMRG results, giving a function A(U, kp) characterizing 

the correction bg^ 2 (L c ) which first increases before de- 
creasing as U varies. When U is small enough, the HF 
approximation reproduces quantitatively the DMRG re- 
sults: Anp(U, kp) Ap,MRG{U,kp). When U becomes 
larger, the decay of Ap,MRG{U,kp) is not quantitatively 
reproduced by the HF approximation. A more suitable de- 
scription could be obtained using a perturbative approach 
adapted to the limit U ^> th, as used in Ref. [53]. For inter- 
mediate U, there is no simple analytical approach, making 
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0.04 
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where the contribution of the conduction band reads: 



Fig. 11. Function A(U,kp) characterizing the interaction in- 
duced correction Sg^ (L c ) to of (|58[) for two interacting 
nano-systems in series (formula lj59| l) at &f = it/2. The exact 
values obtained from the embedding method and the DMRG 
algorithm are shown by circles, the values obtained from the 
HF approximation by a solid line. 

the use of numerical renormalization methods (NRG [22J, 
[24] or DMRG) necessary. 



5 Temperature dependent scale L T of the 
indirect exchange interactions 

We have shown that the indirect exchange interaction be- 
tween two scatterers gives rise to corrections to the value 
of the quantum conductance g 2 F which slowly decays as L c 
increases at zero temperature. We are going to show that 
this effect is suppressed when the length of the wire cou- 
pling the two scatterers exceeds Lt, which is the scale on 
which an electron propagates at the Fermi velocity dur- 
ing a time H/ksT, i.e. the thermal length characteristic 
of free fermions in one dimension. Since our approach is 
essentially valid in the limit of weak values of U, it is 
sufficient to consider the weak interaction limit discussed 
in subsection 3.1 for a temperature T = 0. In this limit, 
the effect was given by the deviation of (Cp +1 c p ) from its 
uniform value sin/cp/V, deviation induced inside the leads 
by an embedded scatterer. To show that this deviation is 
exponentially suppressed above Lt is enough for proving 
that the effect of the indirect exchange interaction upon 
(?2 F vanishes when L c > L^. 

At finite temperature k B T = the Fermi-Dirac 

function f(E, //p) gives the occupation number of the level 
of energy E at a Fermi chemical potential /Up: 



fp(E,n F ) 



1 



(60) 



When one has a single scatterer embedded in an infinite 
perfect lead, the temperature modifies [2U] the value of 
(c\cq) inside the scatterer. Instead of having Eq. ([23]) . one 
now has: 



L 



o 

w dk 



dkfe{E k ^ F )i^»{imo) 

I i=0 

4u cos k sin 2 k 



the contribution of the two bound states becoming 

^bsV) = £ f (E hs: / i P )vc(i)V'bs(o) 

bs 

1-v- 



■ {fp(Ebsl, Mf) — ff3(Ebs2,HF)} ' 
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Fig. 12. v given by Eq. (|62|l as a function of the temperature 
T — ft" 1 for different values of U and fiF = 0. 



The effective hopping term v is given by the implicit 
equation: 



v = l + U(c{co(v)) p , 



(62) 



^cb°(/?,MF) 



Xs°(Amf) 



(61) 



and becomes dependent on the temperature T and on /ip . 

Assuming a half-filled chain (/ip = due to particle- 
hole symmetry), v given by Eq. (|6"2"|) has been calculated 
numerically for weak values U of the interaction strength. 
The T dependence of v is shown in Figs. One can see 
that this dependence remains weak when the temperature 
k B T < 0.1 and U is small. 
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Fig. 13. 



A as a function of p for U = 0.5 and 



Hf = (half-filling). A is the asymptotic value when p — » oo. 

The equation (^SJ giving (Cp.j.jCp) outside the part 
where the electrons interact (p > 1) becomes at a tem- 
perature T: 



2^ 



i=0 



40 
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Fig. 14. Thermal length Lt as a function of the temperature 
T obtained from the exponential decays of (c f p+1 c p )f3 shown in 
Fig. 1131 The points obtained using two values of U = (0.5 and 
1) and T form a single solid curve Lt ~ l/(7rfcsT). 



' ^2 f0( Ehs > MF^bs (p + 1) ifjbs (p) 



(63) 



bs 



dkfp(E k ,n F ) {cosfc - G(A,«)} 



+ {//3(£ ; b s l,MF) ~ fl3{E hs 2,HF)} 



-2p-l 



with the length on which a free fcrmion of speed vf = 
fir 1 dEp / dk-p — 2/h propagates during a time S/(fcgT) 
in one dimension. Ignoring a multiplicative factor l/(27r), 
Lt is the usual thermal length of free fermions in one 
dimension. 



the function G(k,v) being defined in Eq. (|27p. 

By using numerical integration, we have calculated 
(Cp +1 c p )/3 as a function of p for different temperatures T 
and /xp = 0. The data show that (c, 



: jt>+l c P/>/3 



exhibits oscil- 



lations which have a faster decay when the temperature 
increases. If A is the asymptotic value of (c P+1 Cp)f3 when 
p — > oo, we show in Fig . [T5] how the absolute value of 
the difference between (Cp +1 c p )^ and A decays when p 
increases for different values of T. One can see that the 
decay towards its asymptotic value of {c P+1 c p )fj becomes 
exponential: 



c p+i c p) fj A 



cc exp 



P 

Lt 



(64) 



when p is large. 

In Fig. [TJ] the decay length Lt is shown as a function 
of T. One can see that Lt decays when the temperature 
increases as: 



(65) 



Since we have taken = 1 and the lattice spacing s = 
1 in our calculations, one can identify this decay length 



6 Conclusion 

Extending the Landauer formulation of quantum trans- 
port to nano-systems inside which electrons interact, we 
have studied a one dimensional spinless model which is 
simple enough to be analytically solved assuming the 
Hartree-Fock approximation. We have shown that the 
scattering becomes non local when the many-body ef- 
fects inside the scatterer are taken into account. Using 
two identical nano-systems in which interaction gives rise 
to scattering and which are coupled by a non-interacting 
lead of length L c , the Hartree-Fock approximation have al- 
lowed us to map the many-body scatterers onto effective 
one-body scatterers which depend on the other scatterer 
through the indirect exchange interaction via the conduc- 
tion electrons of the coupling lead. The non local charac- 
ter of the scattering is only due to exchange terms in the 
studied model, a positive background charge compensat- 
ing the Hartree contribution. We have shown that the HF 
theory provides a qualitative understanding of the non lo- 
cal effects found in Ref. [8] , quantitatively reproducing the 
exact behaviors in the weak interaction limit. 

Eventually, we point out that we have restricted our 
HF study to a purely one dimensional limit where power- 
ful renormalization methods (DMRG, NRG) are available. 
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Though it has allowed us to compare the HF results to ex- 
act numerical results, this is the worst limit for using the 
HF approximation. To extend the HF study to the many 
channel case, where the electron dynamics will be more 
two or three dimensional does not present particular dif- 
ficulties [7]. Moreover, the HF approximation is believed 
to become more accurate outside one dimension. A HF 
study of this many channel limit is in progress and will be 
published in another work. 

We believe that the interaction induced enhancement 
of the sensitivity of the quantum conductance to the na- 
ture of the attached leads can be relatively easily observed. 
A more detailed study of a possible experiment will be dis- 
cussed in a separate work. 
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8 Appendix 

For two scatterers in series, one uses even and odd stand- 
ing wave functions. The factors ao and a\ for the even and 
odd functions of energy inside the conduction band read: 



-1/2 



«0 



a 



«i 



Oil 



v\ sin fc| 



(V - l)cos \-{L c + 1) 



+ sinfc {v 2 sinfc + (v 2 - 1) sin(fc(L c + 2))} , 
v\ sin fc| 



(^-l)sin( £(L C + 1) 



+ sinfc {v 2 sinfc - (v 2 - l)sin(fc(i c + 2))} . 

The factors bo and b\ for the even and odd functions of 
energy outside the conduction band read: 



bi = 



e *o/2 cosh go(^±3) 



Using the auxiliary functions A±(L C , K) 

smh K 

the normalization factors Aq and A\ for the even and odd 
bound states of energy below the conduction band read: 

-1/2 



-4 



\/2e iJ W2 



2 + 2 coth K 



v 2 A+(L c ,K ) 



A! 



2 + 2 coth Ki 



v 2 A_(L c ,K 1 ) 



, smh 



/<i(L c +3) 
2 
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